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ON  CUMULATIVE  SUM  PROCEDURES  AND  A  STOPPED  WIENER 
PROCESS  FORMULA 

by 

Rasul  A.  Khan 


Abstract 


Let  X-jV  X2,  . ..,  X^  X£,  X^^*  be  independent  random  variables  such 

that  X^,  XJ»\.  X^  -  are  iid  N(0,a^)  and  X^,  X^+^,  ...  are  iid  N(p,c^), 
o  >  0  ,  where  a  is  known  and  k  is  an  unknown  time  index  of  a  possible  change 

V/’ 

in  distribution  .''ror  detecting  changes  in  jf  three  types  of  cumulative  sum  (cusum) 
procedures  are  considered.  The  first  one  is  a  class  of  cusum-type  procedures  -*■— 

s'"  \  ;<1  ^'e" *1  rr— 

such  that  /  E^x  *  +  />  and  (E  Ti<  ^  for  ft  >  0  .  The  second  is  a  modification  of  the 
conventional  cusum  procedure  of  Page  (1954)  which  is  more  efficient.  The  third  is 
a  continuous  version  T  of  the  modified  cusum  procedure  in  terms  of  Wiener  process 
and  its  Laplace  transform  is  found  which  leads  to  the  kr  vn  results  of  Taylor  (1975) 
and  Nadler  and  Robbins  (1971)  .  - - 


AMS  1970  Subject  Classifications:  Primary  62L10;  Secondary  62L99 


Key  Words  and  Phrases:  PetecXton,  cumuJta, tivi  6um  [ca&um] ,  average  Auk  length  (ARLj, 
W+eneA  paocma,  Lap&icz  Pianino  Am. 
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1.  Introduction 


Random  samples  of  size  m  are  taken  at  regular  intervals  from  a  production 

process  and  sample  means  X^,  X2,  . ..  are  computed.  It  is  assumed  that  X^,  X2 ,  .. 

are  independent  random  variables  having  normal  distribution  with  mean  y  and  known 
2 

variance  a  .  The  mean  y  is  said  to  be  in  control  if  y  *  yg(y<yg)  and  out  of 
control  if  y  >  y^  .  There  is  no  loss  of  generality  in  assuming  that  y^  «  0  . 

Thus  X^,  X2,  ....  X^.  X^.^,  ***  are  assumed  to  be  independent  random  variables 

such  that  X^,  . ...  Xj^  are  N(O.o^)  and  X^,  X^.^,  ...  are  iid  N(y,o^), 

y  >  0  ,  where  a  is  known  and  k  is  an  unknown  time  index  of  a  possible  change  in 
distribution.  The  oldest  method  for  detecting  changes  in  y  is  the  Shewart  (1931) 
control  chart.  Motivated  by  Wald's  sequential  probability  ratio  test  and  a  desire 
for  quick  detection  Page  (1954)  defined  the  following  cumulative  sum  procedure. 


Let  Y. 


X  -  r  ,  r  s  0  ,  and  set  S  =  0  ,  S  -  £  Y.  ,  Wn  -  0  , 


Wn  *  max(0,  +  Y  )  .nil.  For  h  >  0  Page's  (1954)  one-sided  cusum  procedure 

is  defined  by  the  stopping  variable 


(1.1) 


t  -  inftnil:  Wn>h}  =  inf{n:  Sr  -  05j§n  S1>h}  , 


and  a  corrective  action  is  taken  at  Wfcih  . 

The  average  run  length  (ARL)  is  defined  to  be  E^t  before  the  corrective  action 

is  taken  while  the  mean  has  remained  at  a  constant  level  y  .  The  rationale  for  this 

(k) 

definition  is  as  follows.  Let  P  denote  the  probability  under  which  X.. ,  ....  X,  .. 

2^0  ^  K— A 

are  iid  N(0,o  )  and  X,  ,  X,  . ,  ...  are  iid  N(y,cr  )  ,  y  >  0  ,  where  a  is  known. 

(°°)  ^  ^  /I) 

Thus  Pn  ■  P_  *  entails  the  model  of  no  change  and  P  =  Pv  '  means  a  change  right 
°  °’W  (k)  M  Vk^ 

from  the  start.  Let  E  denote  the  expectation  under  Pv  '  .  One  would  like  to 

”  fkV 

define  a  detection  stopping  variable  x  such  that  sup  Ev  ' ( (x-k+1) | x>k-l)  is  minimum 

(«)  ksl  ^ 

Subject  to  E^  x  "  *  where  A  is  a  preassigned  positive  constant.  It  turns  out 

that  t  defined  by  (1.1)  has  the  property 


(1.2) 


sup  E^)  ( (t-k+1)  1  t>k-l)  =  E^t  -  E  t 


To  see  this  we  observe  that 


(1.3) 


sup  E^)  ((t-k+1)  |  t>k-l)  £  E^t 

v>i  w  y 


Next,  note  that  t  can  be  written  as 


n 

t  -  inf{n£l:  max  Z  Y. ih)  , 

0£j<n  i=j+l  1 

and  define  the  cusum  procedure  t^  in  terms  of  Y^,  Y^+^,  ...  as 

t,  *  inf{n£k:  S  -  min  S.£h}  , 

n  k-lSj<n  2 

n 

where  S  =  Z  Y.  =  S  (k)  ,  n£k  . 
n  i-k  1  n 

Then  t  -  inf 
k2l 

sup  E(k)  ((t-k+1) I t>k-l)  s  sup  E(k>  t,  -  E(1)  t,  -  E  t  , 
ka  *  kki  “  k  “  1  “ 

which  combined  with  (1.3)  jusitifes  (1.2)  and  the  definition  of  ARL. 

There  is  vast  amount  of  lieterature  on  the  cusum  procedure  (1.1).  The  constant 
r  is  basically  a  design  constant  so  as  to  minimize  ARL^)  at  a  fixed  y^^  £  0 
subject  to  EQt  ^  A  .  This  problem  has  been  treated  by  Ewan  and  Kemp  (1960)  in  the 
normal  case  and  by  Khan  (1978)  for  the  general  family  of  exponential  distributions. 
Unfortunately  Page’s  cusum  procedure  and  the  Shewhart  control  chart  as  well  as  the 
moving  average  procedure  of  Lai  (1974)  have  finite  ARL  when  y  -  0  .  However,  in 
some  problems  a  more  desirable  property  would  be  E^x  =  +  «  while  E  t  <  «®  for 
y>0  ,  which  should  be  as  small  as  possible  or  at  least  fares  well  when  compared  with 
the  conventional  procedures.  A  trivial  Sq  procedure  with  infinite  ARL  when  y  -  0 
is  given  by 

t.  -  inf (nil:  S  ih}  , 
n  In 

where  S-ZY.  ,  Y  -  X,  -  r  ,  r  >  0  . 

n  1-1  1  1  1 

However,  E^t^  -  +  00  for  0  <  y  <  r  so  that  small  changes  cannot  be  detected  by 
t^  .  thus  it  is  deslreable  to  develop  a  cusum-type  detection  procedure  with  the 
above  mentioned  properties. 

A  summary  of  this  paper  is  in  order.  In  Section  2  we  develop  a  class  of  cusum- 
type  detection  procedures  x  such  that  E^x  -  4*»  and  E^x  <  «  for  y  >  0  .  A 
modification  of  the  cusum  procedure  is  given  in  Section  3  and  there  is  numerical 
evidence  that  the  modified  procedure  is  more  efficient  that  the  conventional  cusum 


(t^+k-1)  S  t^  +  k  -  1  ,  and  hence 


procedure.  Finally,  in  Section  4  we  discuss  a  continuous  version  of  the  modified 
cusum  procedure  in  terms  of  Wiener  process  and  obtain  its  Laplace  transform  which 
leads  to  simple  alternative  derivations  of  some  of  the  results  of  Taylor  (1975) 
and  Nadler  and  Robbins  (1971) . 
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2.  A  Cusum-Type  Procedure  t  with  EqX  “  +°° 

We  will  use  the  likelihood  ratio  and  the  mixing  techniques  of  Robbins  (1970) 

to  define  and  study  a  class  of  cusum-type  procedures  x  with  the  property  E^x  =  +® 

and  E  x  <  «  for  y>0  .  Let  X1 ,  X- ,  . . .  be  independent  normal  random  variables 

with  the  mean  y  and  known  variance  o  .  The  mean  y  is  in  control  if  y  S  0  and 

y  >  0  indicates  lack  of  control.  Clearly,  it  is  enough  to  consider  y  -  0  versus 

(k) 

y  >  0  and  assume  that  a  *  1  .  Let  denote  the  probability  under  which 

Xlt  ....  Xk_1  are  iid  N(0,1)  and  Xk»Xk+1»  ...  are  iid  N(y,l)  (y>0)  random 
variables  where  k  is  an  unknown  time  index  for  a  possible  change  in  distribution. 


Obviously,  (A) 


Pq(A)  if  A  €  BCX^-.-.X^j)  and  P 


(k) 


(A) 


Pw(A) 


if  A  e  BCXj^Xj^,  . . .)  ,  where  P^  denotes  N(y,l)  probability  measure.  If  the 
sequence  X^,...,Xn  is  observed,  its  joint  probability  density  function  under  P^^ 
is  given  by 


k,n  0,n 


n 

n  4>(X. ) 
i-1 


if  n<k 


k-1  n 

-  n  *(X  )  n  $(X,-y)  ,  if  n>k  , 

i-1  1  i-k 

where  <J»  (x)  -  (2ir)  exp  (-x2/2)  . 

Any  sensible  procedure  for  detecting  changes  in  y  would  compare  the  likelihood 

that  a  possible  change  has  occurred  at  some  k(l^ksn)  within  the  observed  segment 

(X. , . . .X  )  versus  the  possibility  that  it  will  occur  in  the  future  (k>n)  .  But 
i  n 

this  means  that  such  a  procedure  must  be  based  on  the  ratio 

Zn,k<“>  ■  £k,n  fo!n  '  ♦«!>  ♦«!> 

-  exp(y  S  ,  -  (1/2) y2  (n-k+1)) 


where  S 


n,k 


n 

2  X 
i-k 


It  is  easy  to  define  a  cusum  procedure  based  on  Z  ,  (y)  if  y  is  known  (c.f.  Khan 

(1979b)).  However,  since  y  is  unknown,  a  possible  approach  is  to  use  a  suitable 

mixing  of  Z  ,  (y)  ,  with  respect  to  a  distribution  function  (df)  F(y)  (cf.  Robbins 
n,K 

(1970)).  This  is  exactly  what  leads  to  a  class  of  cusum-type  procedures  given  below. 


Clearly,  ^(y)  is  a  Pg-martingale  relative  to  Fq  ^  -  B(X^,...,Xn>  ,  n£k 
with  EqZq  ^(y)  «  1  .  Now  define 

'n.k  '  I  Zn,kMiZ<»>  ’ 

*  —CO  9 


where  F(y)  is  a  df  on  (-*,«)  . 

Then  ?  ,  is  also  a  Pn-martingale  relative  to  F  ,  ,  n^k  .  Let  b,  be  ar 
n,k  U  n  j  K  k 

increasing  sequence  of  positive  constants  and  define 
(2.1)  t  *  inf{n2l:  C  .  s  b,  for  some  lsksn) 


If  b^  -  b>0  ,  then  T  reduces  to 

x..  -  inf{nal:  max  ?  b)  , 

1  Ikn  M 

a  procedure  studied  by  Poliak  and  Siegmund  (1975).  However,  our  interest  is  in  t 
which  attains  E^C  -  +*  by  a  proper  choice  of  b^  . 

Since  {?  .  ,  F  .  ,  n£k}  is  a  positive  martingale  with  E_c  .  *  1  ,  it 

HjK  n,K  U  QjK 

follows  from  martingale  inequality  that 

(2’2)  P0(max  k  a  bk1)  *  K1  * 

U  n;>k  n,K  K  k 

Define 


Tk  -  inflnkk:  cn>k  2  bk} 
Thus,  one  obtains  from  (2.2)  that 


P0(x<»)  s  Z  PQ(x,  <»)  -  l  P0(max  C  k  *  b  ) 
k-1  u  k  k-1  0  nak  n’k  k 


s  Z  b”1  s  n  <  1  , 

k-1  K 

k  ®  1  l 

by  a  proper  choice  of  b  ,  e.g.,  with  b  -  a  ,  a  >  1  ,  Z  b,-1  -  (a-1)  <  1  . 

k  k-1  k 

Since  Pq(t<®)  S  n  <  1  ,  EQx  -  +»  .  Moreover,  using  an  argument  in  Poliak  and 
Siegmund  (1975)  it  can  easily  be  seen  that  E  x  <  »  for  y>0  . 


■  %TTi 
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We  choose  ■  b^(a)  ■*■<*>  as  ct  -*■  »  such  as  above  and  all  the  asymptotics 
are  as  a  -*■  »  .  To  obtain  an  asymptotic  upper  bound  for  E^x  define 

Tfc  -  inf {n£k:  ?n>k  >  bk(a)}  , 

and  note  that  x  ■  inf(x.  +  k  -  1)  s  i  ,  so  that  E  t  s  E  x.  .  Assuming  the 

k>l  K  1  y  U  l 

existence  of  F'(y)  it  follows  from  a  result  of  Poliak  and  Siegmund  (1975)  that 

21ogb. (a)  ,  - 

(2‘3)  E/l  Z  f2  lo8  V°>  +  log  ( - 2  - }  “  lo8(2ir(F'(y))  )  -l)/yZ  , 

y 

which  is  really  an  asymptotic  upper  bound  for  E^x  . 

Example  1.  Let  F*(y)  =  $(y)  *  (2ir)  exp(-y2/2)  .  Then 

?n,k  "  <n-k+2)  /  exp(Sn>k/2(n-k+2))  ,  Sn#k  =  X±  , 

2  ___ 

and  taking  bk  *  exp(ak/2)  where  ak  =  i/JkTTogcT  ,  a>l  ,  we  have 


P0(Sn,kaan,k  for  some  n*k)  5  P0(,Sn,kUan,k  for  some  n*k) 

S  exp(-a£/2)  <Tfc  , 

where  an  k  ■  /(n-k+2)  (a^+log (n-k+2) )  . 

In  fact  the  first  probability  is  bounded  by(l/2)exp(-a£/2)  .  Moreoever,  the  one 

sided  cusum-type  procedure  x  reduces  to 

x  ■  inf{nkl:  Ska  .  for  some  l^kSn) 
n,k  n,k 

A  two-sided  cusum-type  procedure  can  be  defined  by 

xn  ■  inf (nil:  |S  |  k  a  .  for  some  lsksn) 
u  n,tc  n,K 


It  follows  from  (2.3)  that  an  asymptotic  upper  bound  for  E^x  ig 
EpT~t2  lo8  a  +  log(— ^ )  +  y2-l]/y2  as  a-*» 

y 


(2.4) 


Let  v.  *  infln2k:  S  ,2a  ,  }  ,  and  since  t  =  inf(v,+k-l)  £  v.  , 

K  n,K  n,K  ksl  k  j. 

E  T  S  E  v  .  Also,  since  a  1  is  an  increasing  concave  boundary,  it  follows 

y  y  •  ®  9  * 

from  a  result  of  Robbins  (1970)  that  an  upper  bound  for  E^v^  is  obtained  by 
solving  the  inequality 


(2.5) 


yVl  5  l0g  “  +  l08(EyVl))  Vl  +  IS)'  +  y  » 


M  y  1  v  ®  U  1  $(y)  * 

where  $(y)  is  the  standard  normal  distribution  function.  The  upper  bound  for 

E  v,  is  in  turn  an  upper  bound  for  E  t 
y  1  y 

Example  2.  Let  F’(y)  *  2<J>(y)  ,  y>0 

*  0  ,  ySO 

g 

In  this  case  ,  -  2(n-k+2)“1/2  - )  exp(s?  ,  /2(n-k+2))  , 

n,k  v£*+2  n,lt 

Ic 

and  with  b^  *  a  (a>l)  (2.1)  reduces  to 

2  8n  k 

t  -  inf{n>l:  |S  ,  I  +  2(n-k+2)  log($( - L=-))  2  a  ,  for  some  lsk<n}  , 

n,K  n’K 

where  a  ,  =  (n-k+2)  log(n-k+2)  +  2(n-k+2)  (k  log  a  -  log  2) 
n»«c 

Moreover,  it  follows  from  (2.3)  that 

(2.6)  E  t  <  [2  log  a  +  log(— — )  +  y2  -  1  -  2  log  2]/y2  as  a  +  ® 

y  y 

The  following  Tables  1  and  2  are  based  on  (2.4)  and  (2.6)  respectively  while 
Table  3  is  based  on  (2.5). 

Table  1 
ARL(y)  (Ex.l) 


y 

/2 

/2 

3/2 

2/2 

4 

2 

2 

71.89 

15.95 

3.54 

2.17 

1.59 

88.69 

20.15 

5.09 

2.64 

1.85 

Table  2 
ARL(u) (Ex. 2) 


_ y _ 

Pi 

/2 

'/l- 

3/2 

2/2 

a 

4 

2 

2 

20 

60.80 

13.18 

3.35 

1.87 

1.42 

50 

77.59 

17.38 

4.40 

2.33 

1.68 

Table  3 


Upper  Bound  for  E^x 

(Ex.l) 

V 

a 

Plfl 

PTfl 

/2 

3/572 

2/2 

20 

86 

21 

5 

3 

2 

50 

105 

25 

6 

3 

2 
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3.  A  Modified  Cusum  Procedure 

The  cusum-type  procedure  t  of  the  last  section  has  the  property  EgX  =  -H» 

and  fairly  reduced  ARL(y)  for  y>0  .  However,  E  x  =  0 which  mean*? 

y  y 

that  it  may  take  a  while  to  detect  small  positive  changes  in  y  -  perhaps  the 
price  for  achieving  EgX  =  +»  .  In  contrast.  Page's  (1954)  cusum  procedure  t  , 

defined  by  (1.1)  with  r  =  0,  has  the  property  that  E  t  *  0(h/y)  and 

2  y 
Egt  ~  h  as  h  -*•  »  (ef.  Khan  (1979a)).  But  the  comparison  favors  x  in  that 

Egt  <  «  .  Thus  it  is  desirable  to  somehow  improve  the  cusum  procedure  (1.1)  to 

increase  Egt  so  that  the  modified  version  becomes  more  efficient.  Later  in 

Section  4  we  will  see  that  this  modified  version  manifests  itself  into  simpler 

proofs  of  continuous  cusum  theory  in  terms  of  a  Wiener  process. 

2 

Let  Y^.Yj,...  be  independent  N(y,cr  )  random  variables  with  known  variance 
a ^  .  As  before  the  target  mean  is  y  *  0  and  y  >  0  indicates  that  the  process 
is  out  of  control.  Let  X.^  *  Yj-k,  k^O  ,  i  *  1,2,...  .  Set  WQ  =  0  and  for  b>0 
define 

Wx  =  0  if  X-^-b  ,  W2=0  if  Wj+X^-b 

*  Xx  if  Xj>-b  ,  -  W1+X2  if  W1+X2>-b  , 


From  Wald's  approximations  (cf.  Khan  (1978))  or  Wiener  process  approximation 
(Section  4)  we  have 


(3.1) 

and 


E  N  » 
V 


1 

(y-k) 


[h  -  b 


(1  ~  exp(-2hy) )  , 
(exp(2by)  -  1) 


y*0 


(3.2)  EqN  =  (h2  +  hb)/a2  ,  y=0  , 

2 

where  y  *  (y-k) /a 

Letting  b  ->  0  these  approximations  reduce  to  the  standard  approximations  to 
ARL(y)  for  cusum  procedure  (1.1).  Choose  k  =  0  (the  purpose  of  the  design 
constant  k  has  been  discussed  earlier  in  Section  1) ,  0=1  and  b  =  h  . 
Then  (3.1)  and  (3.2)  reduce  to 


(3.3) 


E  N(h) 

y 


h _ h_  (1  -  exp(-2hy)) 

y  y  (exp(2hy)  -  1) 


y>0 


and 


(3.4)  EQN(h)  =  2h2 

Letting  denote  Page's  (1954)  cusum  procedure  with  k  **  0  and  boundary  h' 

it  follows  that 


EyVh,)  *  th’  -  (1  -  exp(-2yh'))  ]  ,  y>0 


and  Eq^Oi')  =  h'2 

Setting  h'  ■  h/2  it  follows  that 

(3.5)  EyNj. (hi^2)  4  ~ ~  Ihfi  -  (1  -  exp(-2/2hy))]  ,  y>0 

and  E^Oi/I)  =  EQN(h)  =  2h2 


Thus  N(h)  and  N^(h/2)  have  about  (as  h-**>)  the  same  ARL  when  y=0  while  their 
respective  ARL(y)  are  given  by  (3.3)  and  (3.5).  The  following  tables  compare 
the  ARL(y)  for  N  and  .  These  tables  show  that  at  least  for  large  h  the 
modified  procedure  N  is  more  efficient  that  N.  . 


ARL(y) 


V 

) 

.01 

.02 

.] 

y 

h 

ARL(y) 

0 

.5 

1 

2 

2h2-100 

E^Ch/2) 

100 

18.0 

9 

.5 

5.0 

EuN(h) 

100 

14.13 

7 

.07 

3.54 

2 

E„N1 

590 

46.58 

23 

.79 

12.02 

2h  -590 

y  1 

E  N 
y 

590 

34.35 

17 

00 

H 

• 

8.59 

2 

E„Ni 

940 

60.32 

30 

.16 

15.21 

2h  -940 

y  -L 

E  N 
y 

940 

43.36 

21 

.68 

10.84 

A  Stopped  Wiener  Process  Formula 


A  continuous  version  of  the  modified  cusum  procedure  is  now  consiered  and  in 
addition  to  the  approximation  formulas  used  in  Section  3  we  obtain  simple  alternative 
derivations  of  some  of  the  results  of  continuous  cusum  theory  of  Taylor  (1975)  and 
Nadler  and  Robbins  (1971) . 

Let  (W(t)  ,  W(0)*0  ,  tsO}  be  a  Wiener  process  with  a  drift  parameter  y  and 

scale  parameter  a  .  Let  m(t)  -  W(t)  -  min  W(s)  and  M(t)  -  max  W(s)  -  W(t)  . 

CtesiSt  O^sst 

For  h>0  define 

(4.1)  T-  *  inf{t£0:  W(t)  -  min  W(s)  i  h}  *  inf{t^0:  m(t)  h}  , 

0£s£t 

(4.2)  t-  ■  inf{t£0:  max  W(s)  -  W(t)  2  h}  *  inf{t20:  M(t)  £  h}  , 

Z  O^sst 

(4.3)  and  t  ■  min(x^,X2)  *  inf{tsO:  m(t)  S  h  or  M(t)  £  h} 

Here  x^  and  x^  are  the  continuous  versions  of  Page’s  (1954)  one-sided  cusum 
procedures  while  x  is  a  continuous  version  of  a  symmetric  version  of  two-sided 
cusum  procedure.  Taylor  (1975)  obtained  the  Laplace  transform  of  x^  (hence 
that  of  x1  also)  while  Nadler  and  Robbins  (1971)  obtained  the  Laplace  transform 
of  x  .  Their  methods  are  quite  involved  due  to  obvious  intrinisc  difficulties. 
However,  we  consider  a  continuous  version  T  of  the  modified  cusum  procedure  and 
obtain  its  Laplace  transform  which  lead  to  simple  derivations  of  the  Laplace  trans¬ 
forms  of  x^  ,  and  x . .  In  view  of  the  intrinsic  difficulties  the  methods  used 
here  show  the  power  of  renewal  argument  and  the  straegth  of  Wald's  identity. 

The  continuous  version  of  the  modified  cusum  procedure  is  as  follows. 

For  b>0  and  h>0  define 

(4.4)  T  -  inf{t20:  W(t)  S  -  b  or  W(t)  2  h} 

If  T1  terminates  at  the  lower  boundary  -b  ,  the  Wiener  process  starts  from  zero 

A  I 

all  over  again  and  T^  is  repeated.  The  cycle  continues  until  the  upper  boundary 
h  is  attained.  Thus 

T2  -  inf {tsO:  W(t+T1)  +  b  £  -b  or  W(t+Tx)  +  bJh)  , 

...  ,  etc.  Clearly,  T^  ,  T2  ,  • • .  are  lid  random  variables,  and  the  cycles  of 
T^,  ^2  *  ***  are  rePeate<^  until  the  boundary  h  is  hit.  By  abuse  of  notation 
the  cycle  is  terminated  by  the  auxiliary  geometric  stopping  rule 


(4.5) 


N  -  inf{nil:  W(T  )  ;>  h} 
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and  a  corrective  action  is  taken  at  .  Clearly,  the  run  length  is 


Tx  +  T2  + 


Since  T^,  T^,  ...  are  iid  and  N  has  geometric  distribution,  it  follows  from 
Wald's  lemma  that 

(4.6)  E  T  -  E  T,E  N  -  (E  T,)/P  (W(T.)  £  h) 

p  yip  pl  p  l 

First  we  compute  E  T  by  Wald's  identity  for  Wiener  process  and  then  obtain 

y  12  2 

the  Laplace  transform  of  T  .  Let  a(0)  **  p0  +  -j0  o  .  Then 

iexp(0W(t)  -  t  a(0))  ,  F  ■  B(W(s)  ,  sSt)}  is  a  martingale  with  the  property  that 
E^  exp(0W(t)  -  t  a(0))  *  1  .  It  is  well  known  that  Wald's  identity  holds  for 
defined  by  (4.4).  Thus 

(4.7)  E^exp(0W(T1)  -  Tj  a(0))  -  1 

2 

Set  a(0)  »  0  giving  0  ■»  -2y  ,  where  y  =  p/a  ,  and  (4.7)  gives 

EpexP (”2y  W(TX))  -  1  . 

This  identity  and  the  definition  of  T^  give 

(4.8)  P  (W(T  )  S  h)  »  p  »  (exp(2by)  -  l)/(exp(2by)  -  exp(-2hy))  , 

and  P^(W(T^)  £  -b)  -  q  ■  1-p  -  (1  -  exp(-2hy))/(exp(2by)  -  exp(-2hy)) 

When  y  ■  0  it  is  easy  to  see  that 

(4.9)  Pq(W(T1)  a)  «  P(J  -  b/(b+h)  , 

and  P0(W(T1)  S  b)  ■  qQ  «  -pQ  -  h/(b+h)  . 

Since  E  W(Tj)  -  p  E  W(Tj)  and  EqW2^)  -  a2  EQ  ^  ,  (4.6)  ,  (4.8)  and  (4.9)  give 

E  T  »  Ji _ L_  (1-exp  (-2hy)) 

p  p  p  (exp(2by)  -  1)  »  y*° 

-  (h2+bh)/o2  ,  p-0  . 
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Letting  b  -*■  0  one  finds  the  formula  for  E^t^  (hence  E^t2  also)  which  are 
given  by 

E  T  -  -L-  rh  _  (1-exp, (,-?hy ) j_  ,  _1_  f  (1-exp (2hy))  ,  Q 

yTl  v  1  2y  J  »  y2  y  1  2y  J  ’  * 

2  2 

and  Eq^ i  "  ^oT2  =  h  /o 

We  will  now  find  the  Laplace  transform  of  T  .  To  this  end,  we  determine 
the  conditional  and  unconditional  Laplace  transforms  of  .  Set  a(0)  =  s(sSO) 
and  find  the  two  roots  as 

9+  ■  -y+S  and  6_  ■  - (y+6)  , 

where  S  = 

Hence  it  follows  from  (4.7)  that 

(4.10)  E^jexp(0+W(T1)  -  s  -  E^expO.  W(Tj)  -  s  Tj)  -  1  . 

Let  gx  =  g1(s)  -  E^[exp(-s  T1) 1 W (T1)  -  -b]q  , 

and  g2  -  g2(s)  *  E^[exp(-s  T^)  |W(T^)  *  hjp  , 

where  p(q  ■  1-p)  is  defined  by  (4.8). 

Using  the  definition  of  we  find  from  (4.10)  that 

-b0+  .  h0+  -b0"  .  he-  . 

e  g-j^  +  e  g2  -  e  gx  +  e  «2  "  1  * 

and  after  some  algebra  the  solutions  are 


(4.11) 


Hence 


g^  ■  e  sinh(h6)/sinh(  (b+h)6)  , 


and  g2  *  e  ^  sinh(b<$)/sinh((b+h)6) 


(4.12)  <( >q(s)  ■  E^e  sTl  ■  g^  +  g2  *  (e  slnh(h6)  +  e^Y  sinh(b<5))/sinh((b+h)6) 


t 


where  ,  T2 ,  . . . 


are  lid  and  N  has  a 


Now  recall  that  T  -  T,  +  ...  +  T„ 

1  N 

geometric  distribution  given  by 


P(N*n)  -  pq  ”  ,  n-1,2,...  ,  q-l-p  , 

where  p  is  given  by  (4.8)  if  p*0  and  by  (4.9)  if  p“0  . 
From  (4.11)  we  have 


(4.13) 


E^(e“sll|W(T1)  =  -b)  -  gj/q 


e  bY  sinh(h6) 
q  sinh((b+h)6) 


(4.14) 


P  /_“sTl  uj/t  %  «  ■  ehY  sinh(bd) 

V  JWCTj^)  h)  g2/p  p  sinh((b+h)(S) 


Now  it  follows  from  the  definition  of  T  that 


(4.15) 


E  e-sT  -EE  (e“s(Tl+***Tn)|l^n)  P(N-n) 


Using  conditional  independence  we  have 

)Ep(e“s(Tl+**,Tn)|N-n)  -  Ep(e“sTl|N-n)Ey(e“8T2|N=n) . . . Ep (e“sTn | N=n) 
Moreover,  it  fc'.' ows  from  (4.13)  and  (4.14)  and  the  strong  Markov  property  that 


(4.17) 


Ep(e~sTl|N-n)  -  E  (e"811^^)  -  -b) 


e~by  sinh (h6) 
q  sinh((b+h)6) 


(4.18)  Ep(.-Tn|^)  .  Eu(e-SI"IW(Tn)  -  h)  - 

It  follows  from  (4.15),  (4.16),  (4.17)  and  (4.18)  that 


E  e"sT  -  °T.  q~(n'1)e~(n~1)bY(sinh(h6))n~1  ehYsinh(b6)  .  I 

n-1  (sinh((b+h)6)n“1  *  P  sinh((b+h)6)  ** 


ebYsinh  (b<5 )  . 

sinh((b+h)6) 


e~bYsinh(hd) 

sinh((b+h)6) 
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which  is  the  Laplace  transform  of  defined  by  (4.1).  Since  x2  in  (4.2)  is 
representable  as  if  W(t)  is  replaced  by  -W(t)  ,  replacing  y  by  -y  in 
<J»^(s)  we  obtain  the  Laplace  transform  of  x2  as 


a  result  due  to  Taylor  (1975). 

We  now  turn  to  the  problem  of  finding  the  Laplace  transform  of 
T  -  min(x^,x2)  *  inf{t^Os  m(t)^h  or  M(t)^h)  defined  by  (4.3).  If  m(x)£h  or 
M(t)2ih  ,  then  it  is  easy  to  verify  that  M(x)«0  or  m(x)“0  respectively.  Thus 
M(t^)-0  if  is  the  first  to  stop,  and  m(T2)*0  if  t 2  is  the  first  to  stop. 

Using  this  "starting  from  scratch"  property  and  using  the  argument  of  Khan  (1981) 
we  have 

Lemma  1.  P^T1>T2^  "  Et^ETj+ETj)  ,  P(t1<t2)  -  Et2/(ETj+Ex2)  , 

and  Et  ■  (Et^Et2)/(Et1+Et2) 

Substituting  the  expressions  for  Ex^  and  Et2  (given  earlier)  in  Lemma  1  one 
obtains  the  formula  for  Et  found  by  Nadler  and  Robbins  (1971) . 

Let  <Ks)  ■  Ee  ,<K0)«1,  be  the  Laplace  transform  of  x  .  Using  the  "starting 
from  scratch"  property  and  repeating  the  discrete  argument  of  Khan  (1981)  in  the 
continuous  case  we  have  the  identity 
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a  result  due  to  Nadler  and  Robbins  (1971) . 

Thus  the  continuous  version  of  the  modified  cusum  procedure  gives  the 
approximations  to  ARL(u)  used  In  Section  3  and  provides  simple  derivations 
of  the  Laplace  transforms  of  and  t  which  were  obtained  by  Taylor  (1975) 

and  Nadler  and  Robbins  (1971)  by  difficult  methods. 
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